Faraday Instability in a Surface-Frozen Liquid 
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Faraday surface instability measurements of the critical acceleration, a c , and wavenumber, k c , 
for standing surface waves on a tetracosanol (C24H50) melt exhibit abrupt changes at T 3 = 54 °C, 
~4 °C above the bulk freezing temperature. The measured variations of a c and k c vs. temperature 
and driving frequency are accounted for quantitatively by a hydrodynamic model, revealing a change 
from a free-slip surface flow, generic for a free liquid surface (T > T s ), to a surface-pinned, no-slip 
flow, characteristic of a flow near a wetted solid wall (T < T s ). The change at T s is traced to the 
onset of surface freezing, where the steep velocity gradient in the surface-pinned flow significantly 
increases the viscous dissipation near the surface. 

PACS numbers: 68.10.-m, 47.20.Ma, 61.25.Em, 64.70.Dv 



Spatial confinement of a liquid often changes its prop- 
erties markedly. For example, superheating above the 
equilibrium melting temperature [lj and order quench- 
ing upon freezing [2| were observed under confinement 
only. In particular, flow under confinement is impor- 
tant for processes ranging from tribology to protein fold- 
ing to transport through ion channels in cell membranes 
|3j . A transition from a liquid- like to a granular-solid-like 
shear response was observed at nano-scale confinements 
4]. The surface freezing (SF) effect 5], where a solid 
monolayer forms at the surface of a pure normal-alkane 
(C rl H2„+2) melt, provides a unique system for studying 
semi-confined Sow at a solid-liquid interface. The abrupt 
onset of SF at T s allows one to switch on (and off) the 
solid phase by a small temperature variation. Under- 
standing such interfaces would also elucidate the role of 
flow in nucleation and growth processes of crystals from 
melts, which are dominated by such interfaces 0. 

SF occurs in melts of several chain molecules (alkenes, 
alcohols, semi-fluorinated alkanes, diols, CjEj) 7), and at 
both liquid/liquid |8( and liquid/solid |9( interfaces. Re- 
lated surface ordering effects were observed in melts of 
polymers comprising alkyl chains in the backbone or as 
side chai ns Il0| , in liquid alloys [Tl) , and in several liquid 
crystals |l2l . While the structural and thermodynamic 
aspects of SF have been studied in great detail |l3j , the 
influence of surface ordering on macroscopic near-surface 
flows has received little attention to date in this sizable, 
technologically-important class of materials. 

To study this issue, we employ the Faraday instabil- 
ity, which forms standing wave patterns (SWP) at the 
free surface of a vertically- vibrated liquid [14j . By virtue 
of its simplicity, this instability is outstanding among 
pattern-forming systems and a detailed theoretical de- 
scription has been achieved. The SWP formed depend 
sensitively on, and allow a detailed study of, the changes 



in the surface hydrodynamics upon SF. As this study 
of SF demonstrates, the instability can therefore be em- 
ployed to explore physical processes which are difficult to 
access by other, classical means. It also exemplifies the 
more general class of parametric instabilities, which are 
of interest in a broad range of areas including granular 
media, plasmas, nonlinear optics, reaction-diffusion sys- 
tems, and condensed-matter physics 15]. 

The sample used is n-tetracosane, C24H50. With T s = 
54 °C and a bulk freezing temperature Tf, = 50 °C, it 
exhibits the largest temperature range of SF of all pure 
n-alkanes Our circular, 18-cm diameter, sample cell 
was filled with 300 g of C 2 4H 50 (Aldrich, 99.9% pure) to 
a height of h = 1.5 cm above the cell bottom. The cell 
is immersed in a thermostatted water bath yielding good 
temperature stability (±0.05 °C) and minimizing tem- 
perature gradients between the surface and bulk. The 
bath is vibrated sinusoidally in the vertical direction at a 
driving frequency / by an electromechanical shaker. The 
SWP are tracked and recorded using a CCD camera and 
stroboscopic illumination |l6j . 

Two characteristic quantities were measured: the crit- 
ical acceleration, a c , required for destabilizing the flat 
surface to form the SWP, and the critical wavenumber, 
k c , which defines the spatial periodicity of the SWP at 
the onset of the instability. Square SWP (inset of Fig. 
IHa)) were found for all T and / investigated. Fourier 
transforms were used to determine k c from each pattern. 
The T-dependence of k c and a c at fixed / = 150 Hz is 
shown in Fig. ^ Upon cooling from 58 °C to 54 °C, k c 
remains roughly constant. At 54 °C, however, a ~60 m _1 
jump is observed in k c . Below T = 54 °C, k c increases 
linearly with decreasing T. No cooling/heating hysteresis 
is observed. Similar behavior is seen for a c : it is roughly 
constant for 54°C < T < 58°C, at 54 °C it jumps by 
—20%, and it is roughly constant for T < 54 °C. As the 
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FIG. 1: Faraday instability parameters versus temperature 
for a fixed driving frequency / = 150 Hz. Open (closed) sym- 
bols represent data measured on cooling (heating), (a) Criti- 
cal wavenumber, k c , (b) Critical acceleration, a c . Inset to (a): 
Standing wave pattern at the surface at T = 54° C. Solid and 
dashed lines represent the KM-model calculations assuming 
free-slip and no-slip boundary conditions, resp. Arrows mark 
the discontinuities in the Faraday parameters. 



FIG. 2: Faraday instability parameters measured at 
52 °C (triangles; surface-frozen state) and 58 °C (circles; non- 
surface-frozen state) vs. the driving frequency, /: (a) Critical 
wavelength, k c , (b) Critical acceleration, a c . (c) Dispersion 
relation plot, discussed in the text. The solid and dashed lines 
are calculated from the KM-model assuming free-slip and no- 
slip boundary conditions, respectively. 



discontinuities in both quantities occur right at the onset 
temperature T s = 54 °C reported for SF in C24H50, it is 
reasonable to assign these macroscopic effects to the on- 
set of SF, a conclusion supported by the analysis below. 

The measured fixed- T /-dependence of k c and a c , 
shown in Fig. [21 is monotonically increasing, as is typ- 
ical for viscous Newtonian liquids ^t|j both above and 
below T s . Both quantities are larger for the SF phase 
than for the liquid surface phase over the full /-range in- 
vestigated, in agreement with the fixed-/, T-dependent 
measurements shown in Fig. ^ 

To account for these observations, we turn to surface 
hydrodynamics. The validity of the dispersion relation of 
surface waves on inviscid liquids for accurately predict- 
ing the Faraday instability dispersion relation (fc c vs. f) 
of the SWP has been experimentally demonstrated |l8| 
for liquids of viscosities rj > 100ry wa t e r- Thus, it should 
also hold for liquid C24H50, where ?yc 2 4H 5Q ~ 4?7 wator . For 



a liquid layer of depth h and infinite horizontal extent, 
this dispersion relation is iv 2 = (gk c + 7/pfcj? ) tanh(k c h), 
where g — 9.81 ms -2 , 7, p, and u> are the acceleration of 
gravity, surface tension, density, and angular frequency 
of the wave, respectively 0|. The observed SWP in all 
our experiments exhibit lu = 2nf/2, conforming to the 
common subharmonic response of the Faraday instability 
[l7|. Since tanh(fc c /i) w 1 in all our experiments, a sim- 
ple linear dependence of u> 2 — gk c on k^ is expected, the 
slope of which is j/p. This relation is indeed found both 
above and below T s (Fig. E£[c)), indicating that in both 
the frozen and the non-frozen surface phases the simple 
dispersion relation is obeyed. The larger slope found at 
T = 58 °C as compared to that at T — 52 °C indicates 
that 7(T = 58 °C) > j(T = 52 °C), since p does not 
change upon SF. This is consistent with the standard 
static surface tension curve 7 S (T), measured using the 
Wilhelmy plate method, and shown as a solid line in Fig. 
The onset of SF at T s = 54 °C is clearly manifested 
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in this curve by the abrupt change in the slope from a 
negative to a positive value, due to the drop in the sur- 
face entropy upon SF |5(. Comparing this curve with 
the dynamic surface tension 7d(T), extracted from our 
measurements using the dispersion relation above (Fig. 
13 triangles), reveals a reasonable agreement for the non- 
frozen surface phase (T > T s ). In the SF phase (T <T S ), 
however, a ~ 3 mN/m shift of 7d(7 1 ) below 7 S (T) is ob- 
served, although the linear T-dependence and the slope 
remain the same. 

To account for this downshift, we recall that the pres- 
ence of a surfactant monolayer on a liquid surface has 
been demonstrated to change the dispersion relation of 
surface waves and to increase the damping [jflj . Wc 
have therefore employed a model recently introduced by 
Kumar and Matar for the Faraday instability in a liq- 
uid covered by an insoluble surfactant layer (KM-model) 
[2l| . Without an overlayer, the shear stress at the liq- 
uid surface vanishes, yielding a free-slip boundary con- 
dition for liquid flow. For the surfactant- or SF-layer- 
covered surface, the shear stress is finite, and in the limit 
of large Marangoni number (ratio of surface-tension- 
gradient forces to viscous forces) yields a pinned, no-slip 
boundary condition. We use the KM-model with the for- 
mer boundary condition for T > T s , and with the latter 
boundary condition for T < T s . Given the crystalline, 
solid structure of the surface-frozen monolayer, such an 
assumption appears justified. Using the measured static 
surface tension ~f s (T), the dynamic viscosity rj(T) [22| . 
and the aforementioned change in the boundary condi- 
tion at T s , the KM-model yields for fixed / = 150 Hz the 
lines shown in Fig. ^ ( a ) an d (b) . The good agreement 
with the measured values is evident. In particular, the 
change of the hydrodynamic boundary condition at T s 
results not only in a jump of ^20% in a c , but also in a 
similar jump of ^3% in k c , both in good agreement with 
the experimental observation. Moreover, the frequency 
dependent calculations, shown as lines in Fig. [21 demon- 
strate that the KM-model accounts accurately also for 
the observed frequency dependence of k c and a c over the 
full dynamic range investigated, and for the abrupt in- 
crease in these quantities upon SF. 

The increase in a c upon SF reflects an increase in the 
damping of the surface excitations, assignable, in turn, 
to the large velocity gradient near the surface imposed 
by the pinning of the flow at the surface. The increased 
damping leads to a viscous detuning of the inviscid dis- 
persion relation, and thus to a jump in the spatial pe- 
riodicity of the waves. It is interesting to note that the 
change in the boundary condition is equivalent to an in- 
finite surface elasticity: the surface in our model can still 
deform in the vertical direction, but the deformations are 
such that the surface neither contracts nor expands. 

One of the seminal papers on SF in n-alkanes reports 
light-scattering measurements on thermally-excited cap- 
illary waves above and below T s [23^ . The surface elon- 
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FIG. 3: Static surface tension 7 S measured by the Wilhelmy 
plate method (line) and dynamic surface tension 7d (solid 
triangles) as determined from the surface waves' dispersion 
relation for a fixed / = 150 Hz. 



gations due to capillary waves are microscopic, typically 
of the order of nanometers |24|. However, they are still 
governed by the same dispersion relation as the macro- 
scopic, mm-scale surface waves of the Faraday instability 
[25| . It therefore seems puzzling that for capillary waves 
a decreased damping was observed upon SF, rather than 
the increased damping found here. We believe the dif- 
ferent behavior arises due to the different magnitudes 
of the surface elongations in the two cases, and conse- 
quently the different nature of the excited motion. The 
Faraday instability excites macroscopic waves which pro- 
duce real hydrodynamic flow in the proximity of the 
surface. Thus, the no-slip boundary condition in the 
surface-frozen phase induces an increased velocity gradi- 
ent, and consequently viscous dissipation across macro- 
scopic surface-normal distances ensues. In contrast, for 
the nanometer-sized amplitudes of the thermally-excited 
capillary waves, almost no macroscopic hydrodynamic 
flow occurs. The postulated conformational changes in 
the molecular shape at T s , from a flexible, end-distorted 
shape to a rigid extended one \'2('\ , may even enhance the 
molecular mobility, particularly along their long, surface- 
normal, axis. This, in turn, should reduce the damp- 
ing for these microscopic excitations. Indeed, the in- 
crease in surface-normal molecular mobility, which in- 
creases the surface entropy j^, and the molecular con- 
formation changes, which reduce the internal energy [2fij , 
are among the strongest candidates proposed for explain- 
ing the occurrence of the SF effect. 

The SF-induced changes in the boundary condition, 
demonstrated here to alter significantly the macroscopic 
flow characteristics near the surface for one particular 
alkane, should be of importance for all materials exhibit- 
ing surface or interfacial freezing, particularly when the 
flow geometry involves large surface areas. One example 
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is microfluidic devices [28|. Another example is foams, 
where the thinning and bursting rates of the bubbles are 
determined by the viscous drainage flow within the liq- 
uid films constituting the bubble walls |2^. Thus, the 
present results explain also the hydrodynamics underly- 
ing the increased lifetime of alkane bubbles, reported as a 
macroscopic manifestation of SF [3(| . The increase in the 
effective viscous drag on the drainage flow, caused by the 
no-slip boundary condition in the SF phase, reduces the 
drainage rate, and consequently the wall-thinning and 
bursting rates of the bubbles. This leads to the observed 
increase in the bubbles' lifetime upon the onset of SF. 

It has recently been shown that the Faraday instability 
can be employed to explore rheological behavior of bulk 
liquids j^lj, especially when it is difficult to use other 
techniques (e.g., near a liquid- vapor critical point H^)- 
We demonstrated here that this method allows visual- 
izing and examining on a macroscopic scale rheological 
aspects of microscopic structured surfaces, in particu- 
lar the interesting case of a change of the shear-stress 
boundary condition at a liquid's surface from the generic, 
free-slip one to a no-slip one, characteristic of a wetted 
solid wall . We hope that this study will stimulate 
further experiments focusing on how the surface hydro- 
dynamics is affected by microscopic modifications of the 
structure of surfaces or interfaces, e.g., experiments on 
wetting transitions which have been proven to allow for a 
precise control of liquid surface microstructure as a func- 
tion of temperature jsij. This method could also ex- 
tend the wavelength and frequency ranges of established 
semi-microscopic and microscopic techniques like light 
scatteringj23 and surface x-ray photon correlation spec- 
troscopy [35j towards macroscopic hydrodynamic length 
scales. More generally, the present method could be ap- 
plied to other parametric instabilities to provide deeper 
insight into those areas of science where such instabilities 
occur 
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